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摘 要 : 针对 基于 三 维 视觉 指导 的 运动 想象 脑 机 接口 多 通道 宛 余 信息 较 多 、 分 类 准确 率 差 的 问题 ， 提 出 了 一 种 基于 
小 波 包 分 解 (WPD 一 共 空间 滤波 (CSP) 一 自 适应 差分 进化 《ADE) 的 模式 脑 电 信号 特征 提取 与 选择 分 类 方法 。 首 
先 ， 对 采集 的 多 通道 运动 想象 脑 电信 号 进行 WPD 变化 ， 划 分 出 精细 的 子 频 带 ; 然后 ,分 别 将 WPD 变换 后 的 每 个 子 
空间 作为 CSP 的 输入 ， 得 到 对 应 的 特征 向 量 ; 最 后 ， 使 用 ADE 算法 对 特征 向 量 进行 选择 ， 选 择 出 用 于 分 类 的 最 佳 
特征 子 集 。 采 用 WPD-CSP-ADE 模式 进行 特征 提取 与 选择 ， 较 经 典 的 WPD-CSP 方法 在 分 类 正确 率 、 特 征 个 数 方面 
有 着 更 好 的 表现 。 同 时 ， 所 提 算 法 分 类 性 能 明显 优 于 遗传 算法 、 粒 子 群 算法 。 实 验 结果 表明 ，WPD-CSP-ADE 方法 
能 够 有 效 地 提高 分 类 正确 率 ， 同 时 减少 了 用 于 分 类 的 特征 个 数 。 

关键 词 : 脑 机 接口 ; 运动 想象 ; 脑 电 信号 ; 特征 选择 ; 自 适 应 差分 进化 
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Feature selection algorithm for 3D visual guidance multi-channel motor imagery 


Hu Min, Wang Zhiqiang, Huang Hongcheng!, Li Chong 
(School of Communication & Information Engineering, Chongqing University of Posts & Telecommunications, Chongqing 
400065, China) 


Abstract: Concern the problem that multi-channel Motor Imagery (MI) of Brain-Computer Interface (BCI) based on 3D 


visual guidance with more redundancy information and poor classification accuracy, this paper proposed a pattern 
《 二 classification method based on wavelet packet decomposition(WPD)-common Spatial pattern(CSP)-adaptive differential 
evolution(ADE) for feature extraction of electroencephalogram (EEG) .Firstly, this algorithm used WPD to divide the 
multi-channel motion imagery EEG signals into fine sub-bands. Secondly, it used CSP to obtain the eigenvectors 
corresponding to each subspace of WPD transformation. Finally, it selected the feature vectors through the ADE algorithm 
to obtain the best feature subsets for classification. Using WPD-CSP-ADE mode for feature extraction and selection, it had 
better performance in classification accuracy and number of features than the classic WPD-CSP method. At the same time, 
the classification performance of the proposed algorithm was significantly better than the genetic algorithm and particle 
swarm optimization algorithm. The experiments show that the WPD-CSP-ADE method can effectively improve the 
classification accuracy and reduce the number of features used for classification. 

Key words: brain-computer interface (BCI); motor imagery (MD); electroencephalogram (EEG); feature selection; adaptive 
differential evolution (ADE) 


0 引言 间 模 式 (common spatial pattern，CSP) 对 脑 电 特 征 信 号 进行 

x ARIS, SRMF CSP 缺乏 频率 信息 ， 对 特征 信息 提取 不 彻 
脑 机 接口 (BCI) 作为 一 种 新 兴 的 人 机 交互 方式 ， 在 底 ，Yang 等 人 中 利用 小 波 包 分 解 (C wavelet packet 
脑 和 计算 机 等 外 部 电子 设备 之 间 建 立 了 不 依赖 于 外 周 神经 和 decomposition, WPD) 和 CSP 结合 的 方式 进行 特征 提取 。 
肌肉 组 织 的 信息 传输 通道 03。 基 于 运动 想象 的 脑 机 接口 由 于 WPD 印 作为 一 种 快速 时 频 分 析 算 法 , 能 很 好 地 解决 CSP 存在 
其 自发 性 等 特点 应 用 前 景 广阔 中 。 但 是 ， 如 何 有 效 地 获取 与 的 问题 。WPD-CSP 方法 获得 了 相对 较 好 的 方差 ， 有 效 地 提 
识别 脑 电信 号 仍 是 制约 运动 想象 脑 机 接口 发 展 的 主要 障碍 之 取 到 脑 电 的 空间 特征 信息 ， 但 不 能 很 好 地 去 除 见 余 特征 ， 分 


。 在 脑 电信 号 采集 训练 过 程 中 ， 通 过 视觉 诱导 方式 进行 运 类 性 能 表现 仍然 较 差 。 针 对 这 类 问题 ， 相 关 文 献 中 常 使 用 在 
动 想象 ， 保 证 了 脑 电 信号 的 信息 采集 质量 ， 提 高 了 最 终 分 类 ”大 搜索 空间 中 有 着 优异 表现 的 粒子 群 算法 中、 遗传 算法 (0 等 
WIERZ, 同时 , 为 了 更 好 地 消除 被 试 者 差异 以 及 更 加 准 ” ”进化 算法 对 特征 进行 优化 选取 。 通 过 去 除 见 余 特征 ， 然 后 选 


确 地 识别 被 试 者 的 意图 ， 实 验 通常 会 采集 多 路 通道 的 脑 电信 取 最 佳 的 信号 特征 子 集 ， 进 化 算法 在 很 大 程度 上 提升 了 分 类 

号 。 然 而 ， 特 征 复杂 且 维 度 较 高 的 多 通道 信号 会 消耗 大 量 色 器 的 性 能 ， 提 高 了 分 类 的 正确 率 。 然 而 ， 在 对 脑 电 信号 特征 

计算 资源 ， 同 时 多 通道 信号 存在 过 多 的 见 余 信息 ， 会 干扰 到 ”选择 过 程 中 ， 粒 子 群 算法 忽略 了 脑 电 特征 信号 通道 间 的 相关 

有 效 的 脑 电 特征 ， 造 成 分 类 准确 性 的 下 降 。 性 ， 而 遗传 算法 对 于 信号 特征 的 局 部 能 力 搜 索 较 差 ， 且 其 收 
为 提高 脑 电 信号 分 类 的 准确 性 ， 学 术 界 通过 使 用 共同 空 ” SER See IU. 
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针对 多 通道 视觉 指导 运动 想象 脑 机 接口 中 元 余 信 息 较 多 


和 经 


的 进化 算法 在 特征 选择 


过 程 中 存在 的 问题 ， 


了 利用 WPD-CSP 和 自 适 应 差分 进化 〈ADE) 算法 结合 对 脑 
使 用 小 波 包 分 解 方式 和 
取 ， 对 提取 到 的 特征 集 采 用 自 
取 。 经 过 实验 分 析 ， 本 文 所 提 


电信 号 进行 特征 提取 与 选择 。 


间 模 式 对 脑 电信 号 进行 特征 提 
适应 差分 进化 算法 进行 最 优选 


本 文 提出 


E [n] zr 


特征 选择 算法 相对 有 效 地 提高 了 分 类 正确 率 。 


1 ”特征 提取 方法 
1.1 小 波 包 分 解 


WPD 具有 任意 多 尺度 变化 的 特点 ,将 特征 信号 映射 到 一 
组 互相 正 交 的 小 波 基 函数 张 成 的 空间 上 ， 然 后 细 分 信号 为 高 


频 和 低频 两 个 部 分 ， 


加 精细 的 频率 特征 信息 


各 通道 信号 进行 小 波 分 解 。 


i 层 扩 度 第 j 个 小 波 包 空间 。 


EN 
1 
H 2 


r4 
0 
H, 


接着 对 已 分 


段 设 脑 
对 其 进行 2 层 尺度 空间 划分 , 其 过 程 如 图 


02 1, 考虑 到 Haar 小 波 基 实现 简 六 
运算 速度 较 快 以 及 开销 较 小 , 本 文选 


的 低频 信号 实施 再 次 分 解 ， 
同时 对 高 频 信 号 也 进行 再 分 解 处 理 ， 多 次 细 分 ， 


从 而 获取 更 


] Haar 波 作为 小 波 基 对 


有 信号 的 尺度 空间 为 ， 


£ 


H; 


1 所 示 ，Hi 表示 第 


H? 


图 1 


小 波 


包 分 解 示意 图 


Fig. 1 Schematic diagram of wavelet packet decomposition 


1.2 共同 空间 模式 


CSP 作为 一 种 有 效 的 特 行 


通过 构造 


号 的 作用 


体 实 现 步骤 如 下 : 


方差 矩阵 C 和 C; 。 


a) 计 算 两 类 运动 想象 任务 X, 和 Xo 的 归 一 化 的 协 方差 、 协 


XXi 


F 提 取 手 段 已 被 广泛 应 用 。 
个 空间 滤波 器 ， 使 得 两 类 信和 号 
信和 号 方差 最 大 化 ， 同 时 另 一 类 方差 最 小 化 ， 达 到 区 分 两 类 信 
04 DI, 假设 xi 和 xim 分 别 表示 运动 想象 的 两 类 
任务 ， 其 中 : M 表示 采集 点 数 ，N 表示 通道 数 。 则 CSP 的 


它 是 


经 过 滤波 后 ， 一 类 


trace(XX7) 


XX3 


(1) 


M trace(X,X1) 


trace 表示 对 应 矩阵 的 迹 。 
日 根据 归 一 化 的 协 方差 矩阵 C. 与 C 求 得 混合 空间 协 方 


其 中 : T 是 矩阵 转 


运算 符 ; 


差 矩 阵 C , 


对 混合 空间 协 方差 特征 分 解 : 


C=C +C, -UAUT Q) 


其 中 : 4 IREX AERE 


U 为 特征 向 量 


E 阵 。 同 时 ， 将 


FHUETEDGE FR AB PE A 进行 降序 排列 变换 ， 对 应 的 特征 向 量 和 矩阵 
也 调整 变换 。 
c) 构 造 白 化 变换 矩阵 P ， 进 行 白化 : 
已 = MPU7 G) 
S, =PCPr 
t = PC,Pr a) 
再 对 具有 相同 特征 向 量 的 % 和 呵 特 征 值 分 解 : 
s -UAAUI G) 
S, 2U,AUT 


其 中 : A, 4 AER) ERRER FR KERE ; 


U, ,U; JU Je TH OS IRE 
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向 量 矩 阵 。 同 时 可 得 ,=UV,，4+ 罗 =7 ， 其 中 1/7 为 单位 矩阵 ， 
可 以 得 到 当 Si 对 应 的 特征 值 比较 小 时 ，5; 对 应 的 特征 值 则 较 
大 ， 反 之 也 成 立 。 

dd) 选择 特征 向 量 矩 阵 V，( 或 V。) rni gi m FUfUJE m 7i 
进行 重新 组 合 ， 得 到 新 的 矩阵 Bww ， 然 后 结合 白化 矩阵 ， 校 
建 空间 滤波 器 W 为 

W=B'P (6) 
其 中 : mm 为 正 整 数 ， 同 时 m 的 取 值 范围 为 1<m<[N/2] 。 

e) 利 用 空间 滤波 器 W 对 信号 X 进行 滤波 ， 求 得 新 的 时 间 

序列 z 为 


Z-WX (7) 
将 时 间 序 列 zZ 的 第 j 行 记 为 2Z,(j=43%…2m) ， 同 时 定义 对 
应 方差 v 为 特征 向 量 ， 得 到 特征 向 量 f 。 


v; -lg 
E var(Z;) 
FV] = 


H.rH 


中 :var(2)) 表 示 样 本 Z 对 应 的 方差 。 


2 基于 WPD-CSP-ADE 的 脑 电 信和 号 处 理 


自 适 应 查分 进化 

差分 进化 (differential evolution, DE) 算法 又 称 为 “ 差 
分 演化 算法 ”“ 差 异 演化 算法 ” 由 Store 等 人 04 于 1995 年 提 
出 ， 是 一 种 非常 有 效 的 优化 算法 ， 与 遗传 算法 (genetic 
algorithm, GA)、 粒 子 群 优化 算法 (particle swarm optimization, 
PSO) 一 样 ， 差 分 优化 算法 通过 不 断 的 迭代 ， 进 而 保留 适应 
度 较 高 的 个 体 。DE 算法 的 进化 流程 与 GA 非常 类 似 ， 都 包 
括 初始 化 操作 、 选 择 操作 、 变 异 操 作 和 选择 操作 等 ， 但 其 具 
体操 作 有 所 不 同 ，DE 相 比 于 GA， 其 实施 过 程 中 保留 种 群 的 
全 局 搜索 策略 、 基 于 差分 的 变异 操作 以 及 一 对 一 的 竞争 性 策 
略 ， 在 很 大 程度 上 降低 了 遗传 操作 的 复杂 度 0。 


2.1 


差分 进化 算法 作为 一 个 优化 算法 ， 其 最 终 的 目的 是 求解 
最 优 特征 子 集 : 
min  f(x.x,,.,Xp) (10) 


L H 
st xy Sx Sx 


其 中 : DICGRRBEIRIRIBUAMERES xr. 
x; 取 值 范围 的 上 、 下 界 。 
对 式 (11) 初始 化 操作 ， 随 机 产生 式 (12)。 


x, 分 别 表示 是 第 j 个 分 量 


xix <x;i(0)< x7 

(0) $x, (0) xf TR 
izL2,.,NP;j-L2,...D 

xj, (0) = xf; + rand (0.1) (x75 —25,) (12) 


N 


中 : x, 00 是 特征 向 量 集中 第 0 代 的 第 个 特征 向 量子 集 ; 
n (0) 是 特征 向 量 集中 第 0 代 的 第 i 个 特征 向 量子 集 的 第 j 个 
竺 征 向 量 ; NP 则 代表 特征 向 量 集 的 大 小 。 
a) 变 异 操作 。 差 分 策略 即 随机 选取 特征 向 量 集中 两 个 不 
司 的 特征 向 量子 集 ， 将 其 向 量 差 缩放 ， 然 后 与 待 变异 特征 向 
量子 集 进行 向 量 合成 ， 即 
v(g*1)o xa(g) 

-F(x5(g)-x4(g))iznzn*n 
其 中 ，F 表示 缩放 因子 。 
b) 交 叉 操 作 。 对 第 g 代 特 征 向 量 集合 {%(8)} 及 其 变异 的 
间 特 征 向 量子 集 0n Ce 0] 进行 特征 向 量子 集 间 的 交叉 操 


一 


Yr 


(13) 


Hn 
LI 


站 
nt 
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ven emis "PELLI PLE 
x (8),other 阶 小 波 包 分 解 的 子 频带 个 数 为 > 。 小 波 包 分 解 等 效 为 通过 带 
此 中 ， CR 表示 交叉 概率 ， ji 为 [1.2,…,D] 的 随机 整数 。 宽 相同 的 带 通 滤波 器 ， 则 分 解 的 子 频 带宽 度 为 
o) 选 择 操作 。 利 用 贪 禁 算法 进行 差分 进化 得 到 进入 下 tuu T 
特征 向 量 集合 的 子 集 : E 
uns e e (15) 对 于 经 过 小 波 包 分 解 获取 的 各 个 子 频带 脑 电信 号， 其 获 
x (&).other 取 到 一 定 的 脑 电 频 率 特征 信息 。 但 是 由 于 脑 电信 号 来 自 于 多 


在 DE 算法 中 ，CR 和 5 值 选 择 , 决定 了 特征 向 量 集合 个 采集 通道 , 且 模 式 分 类 任务 与 采集 通道 存在 相关 性 , 因此 ， 
样 性 ， 关 平整 个 算法 局 部 或 全 局 搜索 的 快慢 。 基 础 的 DE $E 使 用 CSP 算法 对 WPD 后 的 各 个 子 空间 频带 进行 特征 提取 ， 

法 中 ，CR 和 FF 值 一 般 为 固定 值 ， 当 CR F 的 值 较 小 时 ， 可 可 以 有 效 地 提取 到 脑 电 信号 的 空间 特征 信息 。 但 是 ， 在 CSP 
以 提高 种 群 局 部 搜索 能 力 ， 加 速算 法 收敛 ; 反之 , 当 CR 与 F 特征 提取 中 , 投影 矩阵 参数 m 值 的 大 小 决定 了 所 选 特征 向 量 
的 值 较 大 时 ， 可 以 提高 种 群 全 局 搜索 能 力 ， 但 是 收敛 速率 较 ”个 数 ， 同 时 也 发 现 投 影 矩 阵 前 m 列 特征 向 量 和 后 1 特征 向 
慢 。 为 此 ， 对 于 每 次 迭代 中 CR 和 下 值 采 用 自 适应 调节 ， 当  ” 量 方差 大 小 决定 了 特征 差异 大 小 , 直接 影响 最 终 的 分 类 识别 。 
种 群 内 的 个 体 适应 度 趋 于 一 致 或 者 收敛 于 局 部 最 优 解 时 ， 而 对 于 m 值 的 确定 , 其 与 脑 电信 号 直接 有 关 ， 并 不 存在 固定 
CR 和 F 适当 增 大 ， 增 加 种 群 的 多 样 性 ， 当 种 群 内 的 个 体 适 值 。 通 常 在 脑 电 信号 处 理 中 ,采用 遍历 MM 值 的 全 部 取 值 可 能 ， 


应 度 较 为 分 散 ， CR 和 F 适当 增 大 ， 加 快 收敛 速率 。 对 于 CR 确定 用 于 该 脑 电信 和 号 的 最 优 值 。 但 是 采用 这 种 操作 方法 ， 导 
TU F 的 自 适 应 调节 如 下 : 致 了 训练 过 程 中 计算 量 消耗 较 大 的 问题 ， 尤 其 对 于 脑 电 采集 
在 每 次 变异 操作 中 ， 对 随机 选择 的 三 个 个 体 进行 从 优 到 ”” 通道 较 多 的 情况 。 而 若 采 取经 验 值 方法 确定 mii. 虽然 可 以 
劣 的 排序 ， 得 到 mx， 其 分 别 对 应 的 适应 度 为 fos fus f s 提高 训练 效率 ， 降 低 计 算 消耗 ， 但 是 由 于 个 体 间 差 异 ， 甚 至 
此 时 变异 操作 变 为 同一 个 体 不 同时 间 段 所 采集 的 脑 电 信号 也 存在 明显 特征 差 
v(g 1) x (8) FQ (8)-x,(8)) (6) ” 异 , 可 能 导致 选取 的 m 值 与 其 最 优 值 存 在 较 大 差异 ， 影 响 到 

其 中 : 5 的 取 值 将 根据 每 次 迭代 生成 差分 向 量 的 两 个 个 体 自 最 终 的 分 类 识别 。 
适应 变化 : 综合 考虑 m 值 选取 以 及 WPD-CSP 特征 提取 中 的 元 余 信 
E E) À (17) 息 ， 都 可 能 直接 影响 到 脑 电 信号 模式 分 类 的 准确 性 ， 为 此 提 
f. f 出 对 已 提取 的 特征 再 利用 ADE 进行 特征 选择 ， 具 体 流 程 如 
其 中 :为 上 限 缩放 因子 ;为 下 限 缩放 因子 。 图 2 所 示 。 首 先 ， 使 用 共同 平均 参考 (common average 
而 交叉 操作 中 交叉 概率 CR 变 为 reference,CAR) 与 数字 滤波 器 对 运动 想象 脑 电信 和 号 进行 预 处 
DEM cop ed uos HOUSSE SUDORE 节律 频段 信号 ;接着 , 根据 节律 频 
CR = Je — fin (18) 段 信 号 的 带宽 以 及 系统 处 理 的 精度 与 开销 ， 选 取 合 适 阶 数 的 
CR, aT WPD 进行 子 频 带 划分 ;然后 ， 对 WPD 每 一 个 子 代 空间 进行 
其 中 : RERE n 的 适应 度 ; fu 表示 当前 种 群 中 的 最 优 ， ”CSP 特征 提取 ， 并 保留 全 部 特征 向 量 。 此 时 ， 通 过 保留 CSP 
个 体 适 应 度 ;， fus 表示 当前 种 群 中 的 最 差 个 体 适 应 度 ，7 了 为” ”滤波 后 全 部 特征 的 操作 , 将 m 值 取 值 问题 转换 为 元 余 信 息 处 
当前 种 群 适应 度 平 均值 ， CRs 为 上 限 交 叉 概 率 ; CR': 为 下 限 ” 理 问 题 ,， 有 效 地 解决 CSP 投影 矩阵 中 复杂 的 m. 值 取 值 问题 ， 


交叉 概率 。 MORE i 
LEX EH HU EE , WA y , Wi 村 
ADE 算法 是 在 DE 算法 操作 流程 基础 上 ,使 用 新 的 变异 UTOR AED HE SOUS N , MRARA 4, NARR 
操作 进行 代 葵 ,并 通 过 式 (17)(18) 确 定 每 次 进化 种 群 的 缩放 因 。 人 征 向 量 个 数 为 
子 与 交叉 概率 值 ，ADE 算法 的 具体 流程 如 下 : PY 52N (0) 
DHE ADE 算法 控制 参数 ， 确 定 适 应 度 函数 。 控 制 参数 : _ "— 
E 寺 征 提 HH 全 部 村 Lr E FV,, 
包括 : MRKD NP. FRAMAT Fes KRRIT F Ai 


上 限 交叉 概率 CR, 、 下 限 交叉 概率 CR, 。 个 特征 向 量 存在 更 多 的 频率 、 空 间 等 维度 宛 余 信息 ， 需 要 对 
b 随 机 产生 初始 种 群 。 其 进行 优化 。 为 此 ， 本 文采 用 ADE 算法 对 WPD-CSP 获取 的 
对 初始 种 群 进行 评价 ， 即 计算 初始 种 群 中 每 个 个 体 的 。 ”特征 进行 选择 ， 利 用 ADE 寻 优 机 制 ， 选 择 出 最 佳 特征 子 集 

适应 度 值 。 再 用 于 最 终 的 分 类 识别 。 通 过 ADE 算法 对 特征 向 量 集 合 的 
d) 判 断 是 否 达 到 终止 条 件 或 进化 代数 最 大 。 若 是 ， 则 终 ”处 理 ， 目 的 是 为 了 减少 用 于 分 类 的 特征 向 量 个 数 ， 减 少 特征 

止 进化 ， 将 得 到 最 佳 个 体 作为 最 优 解 输 出 ， 若 否 ， 继 续 。 向 量 的 宛 余 信息 ， 进 而 达到 提高 分 类 准确 性 的 效果 
e) 按 照 式 (16) 进 行 变异 操作 ， 其 缩放 因子 取 值 由 式 (17) 计 另外 ， 考 虑 到 ADE 算法 实现 中 ， 需 要 对 进化 的 每 一 个 

算 所 得 。 个 体 进行 适应 度 评价 ， 以 获取 当前 个 体 的 适应 度 排序 ， 进 而 
人 ) 按 照 式 (14) 进 行 交叉 操作 ， 得 到 中 间 种 群 ， 其 交叉 概率 完成 后 续 的 变异 、 交 叉 以 及 选择 等 操作 。 由 于 ADE 算法 进 

取 值 由 式 (18) 计 算 所 得 。 行 特 征 选择 是 为 了 选择 出 适应 度 高 特征 子 集 用 于 后 续 分 类 ， 
四 按照 式 (15) 进 行 选择 操作 ， 在 原 种 群 和 中 间 种 群 中 选 ”最 终 目 的 是 为 了 提高 分 类 正确 率 ， 所 以 ADE 算法 的 适应 度 

择 个 体 ， 得 到 新 一 代 种 群 。 评价 函数 可 以 选择 分 类 器 分 类 正确 率 来 评价 。 考 虑 到 ADE 
上 进化 代数 gg+1， 转 至 d). 复 法 中 每 一 个 样本 在 迭代 中 适应 度 值 越 小 ， 实 现 进化 的 可 能 

22 基于 WPD-CSP-ADE 的 特征 提取 与 选择 性 就 越 大 , 因此 将 分 类 器 的 分 类 错误 率 作为 适应 度 评价 标准 。 
假设 对 一 通过 带 通 滤波 器 信号 进行 小 波 包 分 解 ， 该 信号 综合 上 述 分 析 , 使 用 WPD-CSP-ADE 算法 对 运动 想象 及 

电 进行 模式 分 类 的 具体 步骤 如 下 ;使 用 最 近邻 对 测试 集 进行 
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识别 ,得 到 识别 率 ,然后 构造 基于 选择 特征 数 的 适应 度 函数 P dd) 根据 节律 频段 信号 带宽 及 系统 处 理 精度 与 开销 ， 选 取 
为 合适 阶 数 的 WPD 对 训练 集 和 测试 集 进行 子 频 带 划 分 ; 
P= gv) Q1) e) 对 WPD 后 的 每 一 个 子 代 空间 的 训练 集 和 测试 集 进 行 

其 中 : 2 是 控制 因子 , 仿真 设置 为 p=0.5 ; N 是 特征 数 ，R 表 ”CSP 特征 提取 ， 并 保留 全 部 特征 向 量 ; 
示 识 别 率 。 £f HH ADE 算法 对 训练 集 的 特征 向 量 进 行 选择 ， 以 分 类 

a) 将 多 通道 脑 电信 号 分 成 训练 集 和 测试 集 ; 器 分 类 错误 率 为 适应 度 函 数 ， 确 定 不 同 特征 个 数 下 的 最 优 的 

b) 对 训练 集 和 测试 集 的 脑 电 信号 用 CAR 进行 滤波 ; 特征 子 集 ; 

c) 根 据 执行 的 运动 想象 任务 ， 确 定 运动 想象 脑 电 信和 号 的 中 使 用 分 类 器 对 训练 集 不 同 特 征 个 数 下 的 最 优 特征 子 集 
节律 特性 ， 得 出 用 于 分 析 的 脑 电 信号 频率 带宽 ， 并 使 用 数字 ” 进行 训练 ， 得 出 该 分 类 器 的 模式 分 类 模型 ， 并 使 用 该 分 类 模 


滤波 器 对 训练 集 和 测试 集 进行 截取 ; 型 对 测试 集 进行 测试 。 


分 类 错误 率 


分 类 正确 性 较 高 
的 特征 子 集 


MAARTE 


图 2 WPD-CSP-ADE 流程 框图 
Fig.2 WPD-CSP-ADE flowchart 


3 ”数据 处 理 视觉 AFSAR o H TERE 被 试 疲劳 而 降低 数据 的 有 效 
性 和 可 靠 性 , 每 间隔 25 组 被 试 者 休息 5 分 钟 , 同时 保证 每 测 
3.1 数据 集 试 100 组 后 再 休息 20 分 钟 。 如 图 4 所 示 , 给 出 了 实验 范式 中 
实验 中 5 名 被 试 者 (年 龄 分 布 为 24~27 岁 ， 且 都 是 右 利 左 、 右 手 运动 想象 的 三 维 动画 训练 场景 。 
F, ,4 男 1 女 ) 首次 参加 实验 ， 实 验 前 均 未 告知 任何 假设 。 提示 音 
实验 脑 电 设备 主要 为 OpenBCI 脑 电信 号 采集 系统 及 具有 64 
通道 的 脑 电 帆 。 依 据 10/20 脑 电 系统 原理 ， 本 次 脑 电 信号 采 | 准备 T amg | —— «m | 
集 电极 主要 有 FC3. FCz. FC4, C5, C3. C1. Cz, C2, CA, EN EM EM S0 17111 
C6. CPI, CP2, CP3, CPz 和 CP4 电极 点 , 参考 电极 为 Fpz， 图 3 实验 范式 
右 耳 乳 突 处 为 接地 电极 ，250HZ 的 采样 频率 。 被 试 戴 好 脑 电 Fig.3 Experimental paradigm 
V8 Jg , 分别 检测 各 个 电极 的 阻抗 , 确保 各 通道 阻抗 低 于 20kQ， TT PEU XE Ws - 
保持 稳定 , 然后 开始 运动 想象 脑 电 信号 采集 。 采集 过 程 中 ， Fa 
被 试 者 根据 屏幕 中 提示 的 运动 场景 进行 运动 想象 。 (a) 左手 运动 想象 的 三 维 动画 视觉 指导 场景 
实验 的 具体 执行 如 图 3 所 示 , 一 个 实验 周期 共 持续 12 s, TX TRZIMIS E ;» 
其 中 前 7 s 为 训练 时 期 。 具 体 包括 : 前 2 秒 静 止 的 三 维 动画 i- ts Fan 
场景 ; 第 2~3 秒 出 现 箭头 并 指示 方向 ， 提 示 进 行 左 手 或 者 右 (b) 右手 运动 想象 的 三 维 动画 视觉 指导 场景 
手 的 运动 想象 任务 ， 在 第 3 秒 时 产生 提示 音 ， 同 时 出 现 运 动 4 左右 手 运动 想象 的 三 维 动画 视觉 指导 场景 
的 三 维 动画 场景 ， 持 续 时 间 为 4 秒 ， 这 期 间 用 户 进行 运动 想 Fig.4 Three dimensional animation visual guidance scene for left and 
; 后 5 秒 为 实验 间隔 。 利 用 头 戴 式 显示 器 作为 呈现 环境 ， right motor imagery 


象 
在 手 推 盒子 动画 为 视觉 指导 场景 下 每 个 被 试 者 分 别 进行 200 ”3.2 脑 电 信号 处 理 
组 的 运动 想象 脑 电 采 集 ， 实 验 过 程 中 随机 出 现 左手 和 右手 的 为 验证 算法 效果 ， 分 别 对 每 个 被 试 者 的 200 组 脑 电 实验 


已 
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数据 进行 划分 ， 训 练 集 和 测试 集 各 100 组 。 — 
首先 对 脑 电 数据 进行 共同 平均 参考 方法 , 然后 , 使 用 FIR f ”实验 结果 与 分 析 

数字 带 通 滤波 器 对 脑 电信 号 进行 滤波 ， 由 于 左右 手 运动 想象 41 SVM 分 类 器 分 类 结果 
所 诱发 的 最 相关 为 Mu (8-13 Hz) 和 Beta (13-30 Hz) 节律 仿真 实验 中 选用 径 向 基 函 数 构造 的 支持 向 量 机 (support 
频段 的 脑 电信 号 , 所 以 通 带 频率 8~30Hz。 对 预 处 理 后 脑 电 数 vector machine, SVM) 分 类 器 ， 利 用 5x5 交叉 验证 方法 动态 
据 进行 小 波 包 分 解 ， 此 时 脑 电信 号 的 下 限 截止 频率 为 .和 上 — 自 适应 的 确定 最 优 的 参数 值 。 对 于 五 名 被 试 者 ， 分 别 使 用 
限 截止 频率 fa 分 别 为 8Hz 和 30Hz， 使 用 二 阶 小 波 包 分 解 ， “GA、PSO、DEB 及 ADE 方法 选择 不 同 的 特征 个 数 作为 特征 子 


数据 集 将 被 分 别 为 [8,13.5]，[13.5,19]，[19,24.5]，[24.5,30]， 集 ， 测 试 其 分 类 正确 率 ， 结 果 如 图 5 所 示 。 
共 四 个 子 频带 。 对 每 个 子 频 带 的 脑 电信 号 进行 共 空间 模式 滤 从 图 5 中 可 以 看 出 ， 五 名 被 试 者 的 分 类 正确 率 随 着 选择 
波 ， 提 取 对 应 子 频带 的 特征 信息 。 其 中 ， 该 数据 集 的 采集 通 。 个 数 的 增加 ， 分 类 正确 率 也 随 之 提高 ， 当 特征 个 数 达 到 一 定 


道 数 为 15， 全 部 子 频 带 经 过 CSP 处 理 后 保留 全 部 特征 向 量 ， 数量 后 ， 分 类 正确 率 趋 于 稳定 。 为 了 进一步 分 析 ， 在 仅 考 虑 
即 获得 60 维特 征 向 量 。 由 于 宛 余 特征 的 存在 , 直接 进行 分 类  ” 被 试 者 差异 的 情况 ， 对 GA, PSO, DE 以 及 ADE 在 不 同 特 
会 影响 到 最 终 的 分 类 正确 率 。 所 以 对 于 待 处 理 的 特征 向 量 集 ， 征 个 数 下 的 分 类 正确 率 进 行 统计 ， 发 现 五 名 被 试 者 (S1~S5) 
接着 使 用 自 适 应 差分 进化 算法 进行 特征 集 优化 ， 选 择 出 部 分 ”在 特征 个 数 分 别 在 12、10、10、13、8 时 ， 达 到 其 最 高 平均 
适合 于 分 类 的 特征 。 具 体 ADE 算法 的 各 项 参数 如 表 1 所 示 。 ”分 类 正确 率 。 由 此 可 以 得 出 ， 在 相同 的 运动 想象 任务 下 ，1 

同时 ， 为 了 对 比 ADE 算法 用 于 特征 选择 的 效果 ,仿真 ” 于 个 体 间 差异 ， 最 佳 的 用 于 分 类 的 特征 子 集 数量 是 不 同 ， 应 
实验 也 对 GA、PSO 以 及 DE 算法 进行 验证 ， 即 GA、PSO、 根据 具体 情况 进行 设 定 。 另 一 方面 ， 为 了 对 比特 征 选择 前 后 


DE 同样 对 基于 WPD-CSP 特征 提取 后 特征 集 进 行 选择 , 三 种 分 类 性 能 以 及 不 同 特征 选择 算法 间 性 能 的 差异 ， 对 五 名 被 试 

算法 的 具体 参数 设置 如 表 2~4 所 示 。 者 在 不 同 特征 个 数 及 特征 选择 算法 下 最 佳 分 类 正确 率 及 特征 

表 1 ADE 算法 的 参数 描述 数 进行 了 统计 , 统计 结果 如 表 5 所 示 ，CSP, 表示 仅 使 用 CSP 

Table 1 Parameter description of ADE algorithm 进行 特征 提取 ， 未 进行 特征 选择 。 从 表 5 可 以 看 出 ， GA. 

参数 名 参数 值 描述 PSO, DE 和 ADE f CSP 情况 ， 平 均 分 类 正确 率 分 别提 高 

NP 200 种 群 规模 了 5.53%，0.25%，3.27%，8.04%， 且 特征 个 数 分 别 平均 下 

G 100 AUC 降 了 65.0095, 70.0096, 82.5090, 63.7590. [Hj&, PSO 算法 

CRH 0:6 上 上 限 交 叉 概率 在 被 试 S2 和 S3 分 类 正确 率 存 在 低 于 CSP, 的 情况 ， 同 时 也 
CRE o 下 限 交 叉 概 率 发 现 其 较 C5P, 算法 平均 分 类 正确 率 无 太 大 差异 。 

susci 通过 上 述 分 析 可 以 看 出 ， 四 类 特征 选择 算法 可 以 有 效 的 

Kos Ae us LIN 对 全 部 特征 集 进行 选择 , 改善 分 类 准确 性 ， 降 低 分 类 器 开销 。 

与 此 同时 ， 使 用 GA、DE 以 及 ADE 特征 选择 算法 也 比 仅 通 


R2 GA 算法 的 参数 描述 


Table 2 Parameter description of GA algorithm 过 改变 CSP 算法 中 投影 矩阵 参数 的 分 类 效果 好 ,其 可 以 有 效 
参数 名 ZH NOR 地 解决 CSP 算法 中 m 值 调 参 问题 。 
PopulationSize 200 种 群 规模 进一步 对 四 种 特征 选择 算法 进行 对 比 ， 发 现 ADE 的 平 
Generations 100 遗传 代数 均 分 类 正确 率 最 高 ， 相 较 于 GA. PSO 以 及 DE, 平均 分 类 下 
PopulationType 二 进 制 编码 编码 方式 确 率 分 别提 高 了 2.38%, 7.77%, 4.62%, HEN ADE 
FitnessFcn 分 类 错误 率 适应 度 函 数 较 GA 和 DE 无 明显 差异 , 比 PSO 增加 2 个 。 整 体 对 比 来 说 ， 
SelectionFcn 轮 盘 赌 算法 选择 算 子 在 SVM 分 类 器 下 ，ADE 分 类 表现 最 好 , GA、DE 次 之 , PSO 
CrossoverFcn 单 点 交叉 算法 交叉 算 子 分 类 表现 相对 最 差 。 
MutationFcn 均匀 变异 算法 变异 算 子 4.2 LDA 分 类 器 分 类 结果 
PIED 2 下 一 代 的 个 体 数 为 了 更 加 精确 地 对 比 四 种 算法 的 差异 ， 减 少 由 分 类 器 性 
BEES 2 tM 能 的 影响 造成 的 实验 差异 ， 又 选用 了 线性 判别 分 析 Clatent 
X3 PSO 算法 的 参数 掺 述 | dirichlet allocation, LDA) 分 类 器 进行 实验 ， 统 计 五 名 被 试 者 
zERC 3 Parameter guum of PSO "PA 在 四 种 算法 下 的 最 佳 特征 数 及 其 分 类 正确 率 ， 如 表 6 所 示 。 
> I 由 表 6 可 以 看 出 ，ADE 的 平均 分 类 正确 率 仍 是 最 高 ， 相 较 于 
m 200 ETE : GA, PSO 以 及 DE， 平 均 分 类 正确 率 分 别提 高 了 2.696. 
MaxDT 100 最 大 迭代 次 数 A n " 
n PEE 惯性 权重 11.52%、3.40%， 特 征 个 数 ADE 较 GA、PS 以 及 DE 分 别 增 
ci 14962 学 习 因 子 1 加 0.8、3、1.2 个 。 
C2 1.4962 学 习 因 子 2 另外 ， 对 比 SVM 和 LDA 作为 分 类 器 ， 四 种 算法 的 性 能 
fitness 分 类 错误 率 适应 度 函 数 差异 。 除 在 DE 算法 下 ，LDA 整体 分 类 正确 率 略 高 于 SVM 
表 4 DE 算法 的 参数 描述 分 类 正确 率 0.2%。 其 他 情况 下 , GA, PSO 以 及 ADE 在 SVM 
Table 4 Parameter description of DE algorithm 的 平均 分 类 正确 率 都 优 于 LDA。 整 体 对 比 来 看 ，SVM 的 平 
参数 名 参数 值 描述 均 分 类 正确 率 比 在 LDA 下 提高 了 2.3596, 选择 的 特征 个 数 无 
NP 200 种 群 规模 明显 差异 。 
G 100 迭代 次 数 综合 SVM 和 LDA 的 分 类 情况 ， 可 以 得 出 ADE 用 于 运 
CR 0.5 交叉 概率 动 想象 脑 电 特征 选择 的 分 类 正确 率 表现 最 好 ， 且 选择 的 特征 
F 0.6 缩放 因子 个 数 无 明显 增多 。 使 用 WPD-CSP-ADE 方法 有 效 的 提高 了 分 


pess 分 类 错误 率 适应 度 函 数 类 正确 率 ， 降 低 了 分 类 器 开销 。 
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EC NEN RS SVM 分 类 器 下 ， 五 名 被 试 的 最 佳 分 类 正确 率 及 特征 数 
aen Table 5 The best classification accuracy and number of features of 
JAAR Mad dg 
ii five subjects under SVM classifier 
3 NMHEEEETTEXULL 村 征 数 /个 
1 CSP,- GA PSO DE ADE CSP,, GA PSO DE ADE 
3 -——À S1 83 87 80 84 89 32 11 8 12 13 
saia S2 70 76 69 71 74 32 13 10 15 14 
IADE] S3 86 89 84 90 94 32 8 7 10 10 
NOH D R a'w S4 79 84 83 83 86 32 13 11 9 13 
S5 80 84 83 83 87 32 11 12 10 8 


(a) 被 试 S1 不 同 特征 个 数 下 的 分 类 正确 率 


表 6 LDA 分 类 器 下 ， 五 名 被 试 的 最 佳 分 类 正确 率 及 特征 数 


156; — in m Table 6 The best classification accuracy and number of features of 
five subjects under LDA classifier 
"| mu 最 佳 分 类 正确 率 /% 村 征 数 /个 
BRNS -SA PSO DE ADE GA PSO DE ADE 
" S1 85 8l 8 9I 10 7 IH B 
F s2 744 6 71 B 8 12 12 12 
^ $3 87 80 91 92 12 9 9 14 
S4 82 79 84 86 135 8 9 
$08 4 ia 人 S5 85 77 82 84 10 6 10 9 
(b) 被 试 S2 不 同 特征 个 数 下 的 分 类 正确 率 5 RA 
针对 运动 想象 脑 电 模式 分 类 中 ， 由 于 多 通道 脑 电 特征 复 
S sAn EN 维度 较 高 、 元 余 信息 量 大 导致 的 分 类 正确 性 的 WE. d 
aw APANE ] 算 消 耗 量 大 的 问题 ， 提 出 了 利用 WPD-CSP-ADE 模式 对 脑 电 
人 | 言 号 进行 特征 提取 与 选择 。 ADE 特征 选择 算法 通过 自 适应 调 
车 缩放 因子 和 交叉 概率 ， 在 全 局 搜索 能 力 和 种 群 多 样 性 中 保 
F 1 持平 衡 。 最 后 ， 通 过 对 比 研究 了 GA. PSO, DE 以 及 ADE 
per 用 于 特征 选择 的 性 能 结果 。 实 验 结果 表明 ，ADE 用 于 运动 想 
EN 象 脑 电 特征 选择 的 分 类 正确 率 表 现 最 好 ， 且 选择 的 特征 个 数 
ê E D R n as 无 明显 增多 ， 其 能 够 有 效 提高 分 类 正确 率 ， 降 低 分 类 器 计算 
销 。 


(c) 被 试 S3 不 同 特征 个 数 下 的 分 类 正确 率 
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